Signatures of band-like tunnelling in granular nanowires 
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We explore the problem of tunneling through disorderd nanowires, comprised of a random dis- 
tribution of metallic grains, by means of a many-body model that captures the essential physics 
of the system. The random configuration of grains gives rise to a smooth band-like set of states, 
which mediates current flow through the nanowire. Analytical and numerical calculations show the 
characteristic signature of this unusual band- like transport to be a quadratic variation of the current 
as a function of the applied voltage (i.e. / ~ V 2 ), a variation that is clearly observed in experimental 
studies of Pt/C composite nanowires. 
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In spite of the vast literature that has appeared over 
more than fifty years, interest in the problem of conduc- 
tion in granular systems remains as high as ever, in large 
part because these systems provide a unique opportunity 
to study transport under conditions where many-body in- 
teractions and disorder both exert a non-perturbative in- 
fluence on carrier conduction. Recent work on this prob- 
lem continues to result in predictions of novel phenom- 
ena, such as single-electron charge soliton propagation in 
one-dimensional arrays of metallic grains^ and coherent 
electron tunneling over large distances in systems with 
large intergranular conductance^* 3 " Experimental stud- 
ies have revealed interesting effects, such as intergranu- 
lar superconductivity* that can lead to global supercon- 
ductivity, or perfect insulating behavior, depending on 
the strength of the intergranular Josephson coupling, and 
spin-dependent tunneling in self-assembled superlattices 
of magnetic nanoparticles£ 

This paper explores the problem of tunneling through a 
disordered system, comprised of a large number of metal- 
lic grains with a random distribution of sizes, and inter- 
grain separations, under conditions where hopping via 
the grains is dominated by single-electron tunneling, and 
the Coulomb repulsion between neighbouring grains is 
strong. A theoretical analysis of this system shows that 
its random configuration smears out any features due to 
tunneling via the individual grains, and gives rise instead 
to a smooth band-like set of states through which the cur- 
rent flow is mediated. Analytical calculations (and nu- 
merical simulations), reveal the characteristic signature 
of this unusual band-like transport to be a quadratic vari- 
ation of the current as a function of the applied voltage 
(i.e. I ~ V 2 ), without any well-defined threshold for the 
onset of conduction. Clear evidence for this behavior is 



found in experimental studies of Pt/C nanowires, which 
show a quadratic variation of their current for the entire 
range of voltage studied. We also investigate the temper- 
ature dependence of the current in these nanowires, and 
show that the observed behavior can be accounted for 
by allowing for a temperature-dependent hopping rate 
between the grains. 

Since the fabrication and basic electrical characteri- 
zation of the nanowires that we study was described in 
detail in Ref. |6|, for further details we refer the reader to 
this publication. The nanowires are formed by the tech- 
nique of electron-beam induced deposition, which yields 
uniform structures with a typical diameter of ~ 50 nm. 
The nanowires are actually comprised of a Pt/C compos- 
ite, however, with Pt crystallites a few nm in size that are 
embedded in an amorphous C matrix. To study the elec- 
trical properties of the nanowires, we deposit Ti/Au elec- 
trodes by electron-beam, and optical, lithography, and 
make two-terminal measurements of the I — V character- 
istics of nanowire segments a few /nm in length, 6 

The motivation for our theoretical studies is provided 
by our recent investigations of the electrical characteris- 
tics of Pt/C composite nanowires, realized by the process 
of Electron-Beam-Induced Deposition (EBID). The fab- 
rication of these nanowires has been described in detail 
in Ref. |f| where we also reported preliminary measure- 
ments of their current-voltage characteristics. While the 
EBID process yields uniform structures several microns 
long, with a typical diameter of 50 nm, as we illustrate 
in the upper inset to Fig. ^ the nanowires are actually 
comprised of a Pt/C composite, with Pt crystallites a few 
nm in size that are embedded in an amorphous C matrix. 
For theoretical modeling of the nanowires, we consider 
a one-dimensional array of quantum wells (QWs) sep- 




FIG. 1: Comparison of experimental and theoretical (lines) 
currents for measurements between contacts 1-2 (rings), 1-4 
(dots), and 3-4 (squares), given to ~ 0.05 eV and To ~ 0.3 eV 
(see text), at T — 4.2 K. 



arated by tunnel barriers that model the insulating C 
matrix (see Fig. [2J ■ The QWs are characterised by their 
widths nf W , i = 1, . . . , N, where N is the total number 
of QWs in the array. Due to the finite distance between 
the grains in the Pt/C wire, the QWs are separated by 
tunnel barriers of the widths nfj B , 1 < i < j < N. In the 
experimental structure, the sizes of the Pt grains and the 
distance between them are random numbers, hence the 
parameters nf W and fiJj B are generated from Normal 
distributions N(^ w , a QW ) and N(n TB ,a TB ), respec- 
tively, where jj9 w / TB and a^ w ^ TB are the corresponding 
expectation value and standard deviation. This is moti- 
vated by the large amount of Pt grains in the nanowire. 
Physically, the one-dimensional array of wells describe a 
single chain of Pt grains in the nano-wire. Since the inter- 
actions between the different chains merely extends the 
set of conductive states in the system, this model pro- 
vides a reasonable account of the transport properties. 

The largest energy scales in the model are the on-site 
and intergranular Coulomb interactions (on-site charg- 
ing energy U ~ 1 eV and intergranular Coulomb inter- 
action Uij > 0.3 eV—). Thus, it is motivated to transfer 
the model into diagonal form by introducing the many- 
body operators \p)(q\,® which describe transitions from 
the state \q) to \p). Diagonalization of the model also 
gives a complete freedom in the choice of parameters 
Ui, Uij,tij, where Uj is the intergranular hopping. Ex- 
perimental work has shown a strong gating effect on the 
I—V characteristics of the nanowire^ suggesting that the 
transport through them is of single electron character. 
Hence it is only necessary to account for the energy sep- 
aration between states like \N',n) and \N'±l,n'), where 
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FIG. 2: Upper panel: One-dimensional array of quantum 
wells separated by tunnel barriers, connected to the left/right 
(L/R) reservoirs. Lower panel: Transition energies of the 
QW array in increasing order constituting the band-like set 
of conductive states. 



N' is the equilibrium ground state number of electrons in 
the nanowire whereas n, n', are state labels. The energy 
difference between the states \N',n) and \N' — 1, n'), de- 
noted by Ai, m / , corresponds to the energy of the removed 
electron, whereas the difference between \N' + l,n) and 
\N', n'), denoted by A 2tm '> and A 3nn > = A lnn , + U, cor- 
responds to the energy of the added electron plus the 
intergranular charging energy. Hence, written in terms 
of the many-body eigenbasis, the model for the Pt/C 
nanowire becomes 
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Pt/C 



E E 
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E*' +m \N' + m,n){N' + to, 



where E^ +m is the energy for the corresponding states. 
In the experimental set-up, the Pt/C nanowire is cou- 
pled to external contact reservoirs, here modelled by free 
electron-like particles. Hence the Hamiltonian for the 
nanowire and the reservoirs is given by 

H = E £ k c L\ + Kpt/C 

kaeL,R 

+ E ( VkC kJ N ' + m,n'){N' +m+ l,n\ +H.c.)(l) 

ka,nn 
m=0,— 1 

where c£ CT (k) creates (annihilates) an electron in the 
left /right (L/R) reservoir at the energy e^, and the hy- 
bridisation Vk between states in the reservoirs and the 
Pt/C nanowire are assumed to be fc-independent. 

In the Hubbard I approximation (HIA)j2ul the Green 
function (GF) for the transition \N', n)(N'±, n\ is given 
in this model by 



G nn '(iuj) = 



P, 



A > — P /T72' 



(2) 



where 7 = 1,2,3, depending on whether an electron is 
removed or added to the nanowire, whereas r = - T L + 
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T R , and T L / R is the coupling to the left/right reservoir. 
The spectral weight P in ri > depending on the bias voltage 
and the temperature, appears due to the non-trivial anti- 
commutation relations between the Hubbard operators. 9 
It should be noticed that the GFs for the conducting 
transitions have to be self-consistently solved for each 
temperature and bias voltage applied to the system. 

By employing the given approximation of the GF to 
the formula for the current through the system derived 
in Refs. llOlllL we end up at the expression 
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where = f(u> — (J-l/r) is the Fermi function at 

the chemical potential Hl/r of the left/right reservoir. 
From this expression for the current, we proceed with an- 
alytical calculations by putting the temperature to zero. 
Thus, the integral in Eq. Q becomes 
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The large number of Pt grains in the nanowire, sug- 
gests to replace the sum over n, nl by an integral. For the 
sake of argument, we assume a rectangular distribution 
of the transition energies A 7ran ' £ [A" 
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, A; nax ], thus 
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This integral is analytically solvable and gives 



h(x) — I arctan x 1 dA 1 



-x arctan a: + — log ( J 



where x — (/x— A)/(PT/2). For |x| < 1, that is, for states 
close to the chemical potential, we may Taylor expand 
this function, giving (letting /jl — » eV, [ir — > 0) 
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Here, the behaviour is dominated by the second term 
since Pr < 1, showing that the current varies as a 
quadratic polynomial of the bias voltage. The function 
in Eq. Q indicates a staircase like current for each con- 
ductive transition. Hence, a smooth character of the 
total current will arise when the transition energies are 
smoothly spread out in the energy interval corresponding 
to the bias voltage applied, with only a few state around 
the equilibrium chemical potential. This suggests that 



the transition energy spacing is > 0.8 eV^ for small 
intergranular hopping. 

In Fig. ^ we plot the measured variation of the current 
in several different nanowire segments (indicated), as a 
function of the square of the bias voltage, and see that the 
resulting data indeed fall closely on a straight line, just 
as suggested by the above analysis. In particular, there 
is no clear evidence for a threshold for conduction, with 
the quadratic variation being observed over the entire 
voltage range. 

While the quadratic variation of the current in Fig. ^is 
strongly suggestive of the band- like transport mechanism 
that we have discussed, we have also performed numer- 
ical calculations of the current from Eq. ©, and these 
provide further support for this mechanism. The lines 
through the experimental data of Fig. ^ represent the 
results of these numerical fits, and agree very well with 
the experimental data. The numerical calculations have 
been performed with just two free parameters (to and 
Tq). The construction of the one-dimensional array of 
QWs separated by tunnel barriers, enables an estimation 
of the hopping parameters = to exp (—fjJj B ), where to 
is a parameter that is equal for all tunnel barriers and for 
all calculations^ The solution of the resulting eigenvalue 
system yields the (transition) energies A 7 „„/ , see Fig. |21 
Finally, the amplitude of the calculated current in Fig. 
^has been fitted to the experimental result by means of 
the coupling T = TqN- 1 / 3 (T l / r = r/2), where T is a 
fixed parameter accounting for the sum of all chains of 
Pt grains in the nano-wire. 

The temperature dependence of the current can to 
good agreement be found partly through the spectral 



weights, P 7 „ n /(T) = 1 - reexp (— iz/f/JA^ 



where 



/3 = ksT, from self-consistent calculations of the GFs 
for finite temperatures, and partly by imposing that the 
hopping parameter goes like to(T) ~ 1 + T/Tq. The 
latter behaviour is reasonable since the tunnel barriers 
between the Pt grains in the C matrix become effectively 
more transparent for higher temperatures. In Fig. [21 we 
display a typical comparison between our experimental 
and theoretical (solid lines) results for the temperature 
dependence of the I — V characteristics. For T > 175 K, 
leakage currents through the substrate layer are expected 
to be significant in the experiment, thus we have omitted 
any calculations for such temperatures. Nevertheless, it 
can be seen in Fig. [3] that the theoretical results are in 
general good agreement with the experimental. 

The use of the many-body states in the present model, 
suggests that we consider the entire Pt/C nanowire as 
an effective interaction region where electrons can be re- 
moved or added via the transitions \N', n')(N' ± 1, n\ at 
the corresponding energies A 7 „ n /. However, the band- 
like transition energies should not be regarded as energy 
bands in the usual sense, where an electron at a given 
energy can propagate freely through the nanowire due to 
the extension of the energy band throughout the region. 
Rather, addition/subtraction of an electron to/from the 
interaction region is followed by a redistribution of the 
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FIG. 3: Theoretical (lines) vs. experimental (symbols) tem- 
perature dependence of the I — V characteristics (measured 
between contacts 1-4), with k ~ 0.75 and v ~ 0.7. 

initial state such that the energy for the transition be- 
tween the initial and final state not exceeds the potential 
drop over the nanowire given by the bias voltage applied. 

In implementing our theoretical model, we have made 
use of a few approximations in order to obtain the excel- 
lent agreement with experiment. The transport proper- 
ties of the conductive states in the nanowire have been 
treated using mean-field theory, i.e. within the HIA. 
This is motivated by the large number of Pt grains 
that contribute to transport through the nanowires. As 
transmission electron microscopy has demonstrated, the 
grains are arranged closely together, although not in ther- 
mal contact, and the hopping of electrons between them 
should tend to smear out the individual character of each 
grain. Hence, possible charge fluctuations of the individ- 



ual grains are negligible. Although the transition energies 
should be renormalized due to strong electron-electron 
interactions in the grains, such a shift can be included 
into the mean-field treatemenfcS and will not cause any 
qualitative change in the nanowire transport properties. 

Treating the temperature dependence of the hopping in 
linear order is motivated by the separation of the energy 
"bands" (> 0.4 eV), which is much larger than the ther- 
mal energy at the temperatures considered here. A more 
elaborate functional dependence on temperature would 
possibly alter our results quantitatively, but would not 
change the general character of the current. 

In conclusion, we have explored the problem of tun- 
nelling through disordered nanowires, comprised of a ran- 
dom distribution of metallic grains, by means of a many- 
body model that captures the essential physics of the 
system; strong on-site and intergranular electron-electron 
interactions. Our analysis shows that the random config- 
uration of the grains smears out any features due to tun- 
nelling via the individual grains, and gives rise instead 
to a smooth band-like set of states through which the 
current flow is mediated. Analytical calculations (and 
numerical simulations, involving only a small number of 
parameters), reveal the characteristic signature of this 
unusual band-like transport to be a quadratic variation 
of the current as a function of the applied voltage (i.e. 
I ~ V 2 ), without any threshold for the onset of con- 
duction. Clear evidence for this behaviour is found in 
experimental studies of Pt/C nanowires, which show a 
quadratic variation of their current for the entire range 
of voltage studied. By imposing a linear dependence of 
the intergranular hopping strength on temperature, we 
have been able to obtain a temperature-dependent vari- 
ation of the nanowire current that is very close to that 
found in experiment. 
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Although the band-like set of transitions depicted in 
Fig. |21 does not represent all the possible transitions 
in the nanowire, these are expected to give the major 
contribution to the conduction. The expansion di a = 
J2 pq {p\dio\q}\p) {q\, where d i<7 annihilates an electron in the 
ith QW, and \p), \q), are eigenstates of the interacting re- 
gion, contains the transition matrix elements (p\di a \q) of 
which a large number have a negligible probability ampli- 
tude. 

The exponential decay justifies a simplification of the 
model to nearest-neighbour interactions in the intergran- 
ular hopping. For better accuracy, however, we have used 



up to the fourth nearest neighbour hopping in our calcu- lations. 



